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1 Introduction 

A classical scheme used to explain energy cascade in turbulence, see e.g. |T2] 
and [3], is based on the picture of the fluid as composed of eddies of various 
sizes. Larger eddies split into smaller ones because of dynamical instabilities and 
transfer their kinetic energy from their scale to the one of the smaller eddies. 
One can think of a tree-like structure where nodes are eddies; any substructure 
father-offsprings, where we denote the father by j £ J (J the set of nodes) and 
the set of offspring by Oj , corresponds to an eddy j and the set Oj of smaller 
eddies produced by j by instability. In the simplest possible picture, eddies 
belong to specified discrete levels, generations: level is made of the largest 
eddy, level 1 of the eddies produced by level zero, and so on. The generation of 
eddy j may be denoted by Denote also the father of eddy j by j. 

Phenomenologically, we associate to any eddy j a non-negative intensity 
Xj(t), at time t, such that the kinetic energy of eddy j is X?(t). We relate 
intensities by a differential rule, which prescribes that the intensity of eddy j 
increases because of a flux of energy from j to j and decreases because of a flux 
of energy from j to its set of offspring Oj . We choose the rule 

keo-j 

where the coefficients Cj are positive. 

This model has been introduced by Katz and Pavlovic [TU] as a simplified 
wavelet description of Euler equations, suitable for understanding the energy 
cascade. The coefficients Cj — 2 Q I J I represent in our model the speed of the 
energy flow from an eddy to its children. The coefficient a is an approximation, 
averaged in time and space, of the rate of this speed. Regarding solutions 
of Euler equations in dimension 3, it may happen (usually as a short term 
phenomenon) that this speed is higher or lower, sometimes that the process 
itself is reversed, that is the energy flows from the smaller eddies to the bigger 
ones: this is known as intermittency. In [5] and [B] it is shown using Bernstein's 
inequality that the rate /3 for the dyadic 3D Euler model lies in the interval [1, |] 
which corresponds to a € [|,4] for the tree dyadic model. As explained in 
section [O] the order of magnitude of Cj that correspond to K41 is: 



The tree dyadic model ([T]) is a more structured version of the so called 
dyadic model of turbulence. The latter is based on variables Y n which represent 
a cumulative intensity of shell n (shell in Fourier or wavelet space) n = 0, 1, 2, . . . 
Here, on the contrary, shell n is described by a set of variables, all X, 's with 
\j\ = n, the different intensities of eddies of generation n. The equations for Y n 
have the form 

Model ([T]) is thus a little bit more realistic than , although it is still extremely 
idealized with respect to the true Fourier description of Euler equations. 

All these models are formally conservative: the global kinetic energy £{t) — 
J2jXj(t), or £(t) — J2 n ^n(t) depending on the case, is formally constant in 
time; it can be easily seen in both cases, using the telescoping structure of the 

series - . However, in previous papers ([§], 0) it has been shown that the 

dyadic model ((3]) is not rigorously conservative: anomalous dissipation occurs. 
The flux of energy to high values of n becomes so fast after some time of evolution 
that, in finite time, part of the energy escapes to infinity in n. 

The same question for the tree dyadic model (TTJ) is more difficult. Intuitively, 
it is not clear what to expect. Even if the global flux from a generation to the 
next one behaves similarly to the shell case ^ , energy may split between eddies 
of the same generation, which increase exponentially in number. Hence there 
is a lot of "space" (a lot of eddies) to accommodate the large amount of energy 
which comes from progenitors in the cascade. 

The main result of this paper, Theorem 12.11 is the proof of anomalous dis- 
sipation also for model (TTJ. To be precise, we have dissipation for a class of 
coefficients Cj which covers ([2]). The proof is similar to the one in [2] but re- 
quires new ideas and ingredients. 

Apart from anomalous dissipation, we consider also stationary solutions, 
showing the existence and uniqueness of such solutions in Theorems l2.2l and l2.3l 
This kind of argument allows and requires a more general model to be studied, 
namely, one needs to insert a forcing term (to find nontrivial stationary solu- 
tions) and we are able to treat also the viscous analogous of the tree dyadic 
model, adding the viscosity term — v2^^Xj to equation (JTJ. The most general 
model that we introduce is thus system ([7}. 

In Section [5] we describe the model and give a short summary of the main 
results of the paper. 

In Section [3] we discuss elementary properties of the model and prove the 
existence of finite energy solutions. 

In Section 2] we exploit the connection between the "classic" dyadic model 
on naturals and the tree dyadic model. If the number of children is constant for 
every node in the tree, then from each solution on the former one can build a 
"lifted" version on the tree which is a solution of the latter. 

Section [5] is devoted to the proof of the anomalous dissipation Theorem 12.11 
in the inviscid unforced case. Self-similar solutions are also discussed. 

In Section [5] we study the stationary solutions. We prove existence and 
uniqueness of stationary solutions of classic and tree forced systems ([SJ and ([7]) 
with and without viscosity. Here the positive force / is required because other- 
wise the unique non-negative stationary solution is the null one. 
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1.1 The decay of Xj corresponding to K41 and anomalous 
dissipation 

In the case of the classic dyadic model ^ , Kolmogorov inertial range spectrum 
reads 

Y ~ k- 1 / 3 

J-n n-jj 

The exponent is intuitive in such case. For the tree dyadic model ([1} the correct 
exponent may look unfamiliar and thus we give a heuristic derivation of it. The 
result is that Kolmogorov inertial range spectrum corresponds to 

(4) 

K41 theory [T2] states that, if u (x) is the velocity of the turbulent fluid at 
position x and the expected value E is suitably understood (for instance if we 
analyze a time- stationary regime), one has 

E[\u(x)-u(y)\ 2 ] ~\x-y\ 2 ' s 

when x and y are very close each other (but not too close). Very vaguely this 
means 

MaO-uMMs-tfl 1 / 3 . 

Following Katz-Pavlovic [TU], let us think that u(x) may be written in a ba- 
sis (wj) (which are norm-one vectors in L 2 ) as 

u(x) = ^X j w j {x). 

3 

The vector field Wj (x) corresponds to the velocity field of eddy j. Let us assume 
that eddy j has a support Qj of the order of a cube of side 2~ Given j, 
take x,y G Qj. When we compute u(x) — u(y) we use the approximation 
u(x) = XjWj{x), u(y) — XjWj(y). Then 

\u(x) ~ u{y)\ = \Xj\ \wj(x) ~ Wj{y)\ 

namely 

\X 3 \\w 3 (x)-w 3 (y)\ ~ Ix-yl 1 ' 3 , x,y G Qj. 

We consider reasonably correct this approximation when x, y G Qj have a dis- 
tance of the order of 2 _ '- , l, otherwise we should use smaller eddies in this ap- 
proximation. Thus we have 

\X j \\w j (x)-w j (y)\~2-M, x,veQj,\x-y\~2-W. (5) 

Moreover, we have 

\w j (x)-w j (y)\=\Vw j (Z)\\x-y\ (6) 

for some point £ between x and y (to be precise, the mean value theorem must 
be applied to each component of the vector valued function wj). Recall that 
J Wj(x) 2 dx = 1, hence the typical size Sj of Wj in Qj can be guessed from 
s 22-3|j| ^ ^ namely Sj ~ 2%^\ , Since Wj has variations of order Sj at distance 
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2 , we deduce that the typical values of Vwj in Qj have the order 2 2 1 I / 2 I- 7 ' = 
2^1. Thus, from ©, 

\wj(x) -Wj(y)\ ~ 2*1^2-^1. 
Along with |[5j this gives us 

1^12*^2-^1 - 2"^ 

namely 

|X,-|~2<-* + i-*>lil =2-tUI. 

We have established (0}, on a heuristic ground of course. 

Let us give a heuristic explanation of the fact that, when anomalous dissi- 
pation occurs, the decay Q appears. In a sense, this may be seen as a confir- 
mation that (j4]) is the correct decay corresponding to K41. Let us start from 
equations (fTJ with Cj ~ 25 Wl, the Katz-Pavlovic prescription. Let £ n be the 
energy up to generation n: 

lil<« 

Then, as will be seen later with equation ©, 

f L = - 2l(n+1) E 

fc|=n+l 

In order to have anomalous dissipation, we should have 

dt r 

If we assume a power decay 

X,-~2"' 7lil . 



Then, since the cardinality of {Qj : \j\ — n} should be of the order of 2 



3?> 



2§(n+i) A|Xfc ~ 25™2 3 "2 _3,? ™ = 2^ _3 ' ? ^™ 

|*|=n+l 



and thus T) = ^r. 



2 Model and main results 

Let J be the set of nodes. Inside J we identify one special node, called root 
or ancestor of the tree, which is denoted by 0. For all j £ J we define the 
generation number \j\ £ N (such that |0| = 0), the set of offsprings of j, denoted 
by Oj C J, such that \k\ = \j\ + 1 for all k £ Oj and a unique parent J with 
j £ Oj. The root has no parent inside J, but with slight notation abuse we 
will nevertheless use the symbol when needed. 

For sake of simplicity we will suppose throughout the paper that the cardi- 
nality of Oj is constant, JjCj =: iV* for all j £ J, but some results can be easily 
generalized at least to the case where %0$ is positive and uniformly bounded. 
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It will turn out to be very important to compare iV* to some coefficients of 
the model. To this end we set also a :— h log 2 iV» so that — 2 2a 

The dynamics of the tree dyadic model is described by a family (Xj)j^j of 
functions Xj : [0, oo) — > R. Its general formulation is described by the equations 
below. (Notice that X$ does not belong to the family and merely represents a 
convenient symbolic alias for the constant forcing term.) 

Xo(t) = f 

jX 3 = -udjXj + Cj X? - CkXjXk, Vj e J (7) 

Here we suppose that / > 0, v > 0, and that the other coefficients have an 
exponential behavior, namely Cj — 2 a l J *l, dj = 2" |r l J 'l with a > and 7 > 0. 

If / = we call the system unforced, if v = we call it inviscid. 

This system will usually come with an initial condition which will be denoted 
by X° — (Xj)j£j. One natural space for X(t) to live is Z 2 (J;R), which we will 
simply denote by I 2 , the setting being understood. The I norm will be simply 
denoted by | • 1 1 . 

Definition 1. Given X° £ R J , we call componentwise solution of system ([7]) 
with initial condition X° any family X — [Xj)j^j of continuously differentiable 
functions Xj : [0, 00) — > R such that X(0) = X° and all equations in system ([7]) 
are satisfied. If moreover X(t) € I 2 for all t > 0, we call it an I 2 solution. 
We say that a solution is positive if Xj(t) > for all j E J and t > 0. 

Existence of positive I 2 solutions is classical and can be found in Section O 
while uniqueness is an open problem. 

This system of equations is locally conservative, in the sense made rigorous 
by Proposition ^ . 2l below. where the following energy balance inequality is proven 

\\X(t)\\ 2 <\\X{s)\\ 2 + 2f 2 f X (u)du-2isY,dj [ X 2 (u)du 

jeJ -> s 

It turns out that in some cases this is in fact an equality and in some cases it is 
a strict inequality. When the latter happens we say that anomalous dissipation 
occurs. 

The main results of the paper deal with anomalous dissipation and stationary 
solutions. 

Anomalous dissipation on the inviscid, unforced tree dyadic model. 

The proof of the next result is given in Section [SJ 

Theorem 2.1. Let %Oj = 2 2a for all j. Suppose a < a and f = v = in 

equations (UJ. Let X be any positive I 2 solution with initial condition X . Then 
there exists C > 0, depending only on \\X \\, such that for all t > 

£(t) :=\\X(t)\\ 2 :=J2^(t)<^ . 
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This theorem holds also if we use the weaker hypothesis 1 < \Oj < 2 2a 
for all j. The statement tells us that the energy of the system goes to zero at 
least as fast as t~ 2 . In Section IQl we show that for this model there are some 
self-similar solutions and that their energy goes to zero exactly like t~ 2 . So the 
estimate of Theorem 12.11 cannot be improved much. 

Stationary solutions for the forced classic dyadic model. 

It will be important for our purposes to switch between the tree dyadic model 
and the classic one, where J is simply the set of non-negative integers with 
Oj :={j + l}for allj. 

To avoid confusion we will use different symbols for the classic system, whose 
equations are the following. 

Y- X {t) = f 

d (8) 
— Y n = -vl n Y n + k n Y 2 _ x - k n+ {Y n Y n+ i, Vrt > 

with / > 0, v > 0, k n = 2^ n , l n = 2T", /3 > and 7 > 0. 

When this model is interpreted as a special case of ((?]) we will have iV* = 1, 
(5 = and /3 = a. Observe that the definitions of solutions given on the tree 
model extend easily to this one, but notice that in this setting l 2 will correspond 
to the standard space of sequences. 

The following theorem deals with stationary solutions, namely solutions con- 
stant in time. We do not detail the proof, since, by what we said above, it is a 
special case of the analogous statement for the tree dyadic model, Theorem [23] 
which is proven in Section [5] 

Theorem 2.2. If f > 0, then there exists a unique I 2 positive solution Y of 
system © which is stationary. Moreover 

ifv = Q then Y n (t) := / • 2-^ n+1 \- 

if v > and 3j > 2/3, the stationary solution is conservative and regular, in 
that for all real s, E„[ 2S " r «(t)] 2 < oo; 

if V > and 3"f < 2/3, there exists C > such that for all f > C the invariant 
solution of ([8]) is not regular and exhibits anomalous dissipation. 

In the inviscid case, this theorem extends an analogue result of [6] where it is 
proved for fi = ~. In the viscous case it extends a result of [S], in which existence 
and uniqueness of stationary solutions are proved for 7 = 2 and /? € f|, | ] . 

Stationary solutions for the forced tree dyadic model. 

An analogous of Theorem holds for the tree dyadic model too. This is proved 
in Section [5] 

Theorem 2.3. Let §Oj = 2 2a for all j. Suppose a < a and f > in equa- 
tions ([7]). Then there exists a unique I 2 positive solution X which is stationary. 
Moreover 

ifv = then Xj(t) := f ■ 2" JiJ 3 ti ( 2a + a! ) for all j 6 J; 
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if v > and < a — a < ^-y, the stationary solution is conservative and 
regular, in that for all real s, Y^jeA^ Xj{t)] 2 < oo; 

if v > and a~a > §7, there exists C > such that for all f > C the invariant 
solution of |7|) is not regular and exhibits anomalous dissipation. 

3 Elementary properties 

We will provide, in this section, some basic results on the tree dyadic model. 
The results are analogous to those provided for the dyadic model in [2] and [5], 
but the proofs require some new ideas to cope with the more general structure. 

We will suppose throughout the paper that the initial condition X° is in I 2 
and that X® > for all j G J. It will turn out that this two properties hold 
then for all times. 

Definition 2. For n > — 1, we denote by £ n (t) the total energy on nodes j with 
< n at time t and £(t) the energy of all nodes at time t (which is possibly 
infinite) : 

£ n (t):= £x?(f), £(t):=J2^(t). 

\j\<n jEJ 

Note in particular that £_i = 0. 

We will use very often the derivative of £„, for n > 0, 

p n (t)=2Y / X j ±X jit) 

b'l<« 

= -2v E d s x l + 2 E cjXjXj - 2 E E c ^ x l^ 

b'l<™ lil<" bl<«feeOj 

= -2// E d i x j + 2c X 2 X Q - 2 E c kX~ k Xk 

\j\<n |fc|=n+l 

so we get for all n > 

±£ n (t) = -2v E diXft) + 2fX (t) 2 E c fc X|(t)X fe (t). (9) 

|j|<7t. |fe|=n+l 

Proposition 3.1. // X® > /or aZ/ j, t/ien any componentwise solution is 
positive. If X is in I 2 , any positive componentwise solution is a positive I 2 
solution, in particular for all t > 0, 

£{t) < {£{0) + l)e 2f2t . (10) 
Proof. From the definition of componentwise solution we get that for all j G J 

Xj(t)= X°e-^o^+^'« CfcXfc ( r » dr + / c j X?(s)e- I> d i+^ c " x ^ r )) dr d a (11) 

Jo 

yielding Xj(t) > for all t > and all j G J. 
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Now we turn to the estimates of £(t). In ©, since X k (t) > we have two 
negative contribution which we drop and we use the bound X (t) <X$(t) + l< 
£ n (t) + 1 to get that for all n > 0, 

J t £n{t) < 2f (£ n (t) + 1) 

so by Gronwall lemma £ n {t) + 1 < (£ n {0) + l)e 2 ^ 2 *. Letting n — > oo we ob- 
tain mni). □ 

Proposition 3.2. For any positive I 2 solution X , the following energy balance 
•principle holds, for all < s < t. 

£{t) = 8(a) + 2f [ X (u)du-2v^d j f X 2 (u)du 

-2 1im / V c k Xl(u)X k (u)du (12) 
Js |fe|=n 

where the limit always exists and is non-negative. In particular, for the unforced, 
inviscid (f = v = 0) tree dyadic model, £ is non-increasing. 

Proof. Let < s < t, then by © for all n > 0, 

£ n (t) = £ n (s) -2vJ2 d J [ Xf{u)du + 2/ 2 f X (u)du 

I _• i J S J S 



2 1 H c k Xl(u)X k (u)di 



\k\=n+l 



As n — > oo, since the solution is in I 2 , £ n (s) t £{s) < oo and the same holds for 
t. The viscosity term is a non-decreasing sequence bounded by 

2v d J I x j( u ) du < £ ( s ) + 2 / 2 / X (u)du < oo 

I •! ■ J 8 J S 

so it converges too. Then the border term converges being the sum of converging 
sequences. □ 

Definition 3. We say that a positive I 2 solution X is conservative in [s, t] if 
the limit in (fT!?)) is equal to zero that is if 

£ (t) = £(s) + 2/ 2 / X (u)du - 2v^d 3 I X 2 (u)du 

Otherwise we say that X has anomalous dissipation in [s,t]. 

Theorem 3.3. Let X° G I 2 with X° > for all j G J. Then there exists at 
least a positive I 2 solution with initial condition X° . 



Proof. The proof by finite dimensional approximates is completely classic. Fix 
N > 1 and consider the finite dimensional system 



f t X } = -pdjXj + cjjq - E k eo 3 CkXjX k jeJ,0<\j\<N 

x k (t) = o k e J, \k\ = N + i 

{x s (p) = x$ jeJ,o<\j\<N 



(13) 



for all t > 0. Notice that proposition 13. II is true also for this truncated system 
(with unchanged proof), so there is a unique global solution. (Local existence 
and uniqueness follow from the local Lipschitz continuity of the vector field 
and global existence comes from the bound in (|10p.) We'll denote such unique 
solution by X N . 

Now fix j G J and consider on a bounded interval [0, T] the family (Xj) N>\j\ ■ 
By (|10p we have a strong bound that does not depend on t and N 

|Xf < (£(0) + 1)M T/2 ViV > 1 V* G [0,T] , 

thus the family (Xj , ) N> y\ is uniformly bounded, and by applying the same 
bound to (TT3"]) . equicontinuous. From Arzela-Ascoli theorem, for every j G J 
there exists a sequence (Nj,k)k>i such that [X- 3,h )k converges uniformly to a 
continuous function Xj . By a diagonal procedure we can modify the extraction 
procedure and get a single sequence (Nk)k>i such that for all j G J, X^ k Xj 
uniformly. Now we can pass to the limit as k — > oo in the equation 



Xf k = X° 



» + f [-vdjX* (r) + a (Xf (r)) 2 - Y c i X ? h (r)*? k W 



dr 



and prove that the functions Xj are continuously differentiable and satisfy sys- 
tem ([7]) with initial condition Xj. Continuation from an arbitrary bounded 
time interval to all t > is obvious. Finally, X is a positive I 2 solution by 
Proposition 13.11 □ 

We conclude the section on elementary results by collecting a useful estimate 
on the energy transfer and a statement clarifying that all components are strictly 
positive for t > 0. 

Proposition 3.4. The following properties hold: 
1. If f = 0, for all n>-l 



\k\=n+l 



c k Xl(s)X k (s)ds < £ n (0) (14) 



2. if Xj > for all j s.t. \j\ = M for some M > 0, then Xj(t) > for every 
j s.t. \j\ > M and all t > 0. 

Proof. 1. If n = —1 the inequality is trivially true. If n > 0, by integrating 
equation © with / = 0, we find that 



£n{t) 



1v [ V djX?(s)ds = £ n (0)-2 [ V c k Xl(s)X k ( s )ds 

Jo \j\<n J ° |fe|=n+l 
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The left hand side is non-negative for all t, so taking the limit for t — > oo 
in the right hand side completes the proof. 

2. For \j\ = M we have from 

Xj(t) > X°e~ Io^ d i+^k c k x k{r))dr > q 

Now suppose that for some j E J \ {0}, Xj(t) > for every t > 0. Then 
again by ((TTJ) , 

Xj(t) > f cjX*(a)e- S> d ^^ ^(r ))d r ds > Q 



By induction on |j| > M we have our thesis. □ 

4 Relationship with classic dyadic model 

Recall the differential equations for the tree and classic dyadic models. 
X- (t) = f 

jXj = -vdjXj + CjX? - CkX 3 X k , Vj e J (15) 

keOj 

Y-t(t) = f 

d (!6) 
— Y n = -ul n Y n + - k n+1 Y n Y n+1 , Vn > 

where / > 0, v > and for all n G N and j E J, 

r . — Q b'l — O0n rf. — 07\j\ 1 _ 07™ 

Again we assume that \Oj = iV* = 2 2 " for all j € J, but we stress that for this 
section this is a fundamental hypothesis and not a technical one. 

The following proposition shows that examples of solutions of the tree dyadic 
model (|15p can be obtained by lifting the solutions of the classic dyadic model (|TB|) . 

Proposition 4.1. If Y is a componentwise (resp. I 2 ) solution of (|16p . then 
Xj(t) :— 2~^^ +2 ^ a Yy\ (t) is a componentwise (resp. I 2 ) solution of (|15p with 
a — /3 + a. IfYis positive, so is X . 

Proof. A direct computation shows that X is a componentwise solution. Then 
observe that, for any n > 0, 



J 3 n n 



\J\- 

SO 



|i|<« k<7i 

Positivity is obvious. □ 
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Remark 1. If we consider a fixed, since /3 = a — a, for small values of 
we'll have larger values of (3, and the other way around. That is to say, the less 
offspring every node has, the faster the dynamics will be. 

Remark 2. Let us stress that j3 > when < 2 2a . Since the behavior of the 
solutions of (fT()]) is strongly related to the sign of /?, then the behavior of the 
solutions of (|15[) is strongly connected to the sign of a — a. For example, in the 
classic dyadic there is anomalous dissipation if and only if j3 > 0, and hence in 
the tree dyadic there will be lifted solutions with anomalous dissipation when 
a > a and lifted solutions which are conservative when a < a. 

5 Anomalous dissipation and self-similar solutions 
in the inviscid and unforced case. 

Throughout this section we'll consider system (O in its unforced (/ = 0) and 
inviscid (y — 0) version. 

X- (t) = 

jX = Cj X 2 - J2 CkXjX k , Vj g J (17) 
keOj 

Equation ([9]), that is the derivative of energy up to the n-th generation 
becomes 

j t £ n {t) = -2 J2 c k X 2 k (t)X k (t), n>0 

\k\=n+l 

Since only the border term survives, one would expect it to vanish in the limit 
n — » oo. This can be rigorously proven only if the solution lives in a sufficiently 
regular space, that is to say that Xj goes fast to zero as \j\ — > oo. For the 
classic dyadic Kiselev and Zlatos [TT| proved that solutions that are regular in 
the beginning, stay regular for some time but then lose regularity in finite time. 
Thus our analysis is not restricted to regular solutions, and in fact we will prove 
in this section that for sufficiently large times all solutions dissipate energy. 

Let us give some definitions. Let us denote by jj the energy at time 
in the subtree Tj rooted in j plus all the energy flowing in j from the upper 
generations, 



Let < s < t and define for all j G J 



/>oo 

7i := E + / ^X^ds 

7.,-Tn JO 



m i - |nf Xj (r) 

Lemma 5.1. Let X be a positive I 2 solution of system (jTTJ) . The following 
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inequalities hold for all n > 0. 

£„(*)- £ n _i(s) < E m 2 j <£(0) 

b'l=n 

E 7j < *(0) 

b'l=n 

E X ^ r ) 2 ^ Ti. Vr > 

Proof. The upper bound is obvious, since 

E E X,(«) 2 <£„(«)<£(<)) 

b'l=« b'l=« 

where we used Proposition 13.21 Now let j E J. From (jTTJ) we have for the 
differential of Xj 

keOj 

Let r G [s,t] and integrate on [s,r], yielding 

X?(r) = X?( S ) + [ 2c ] X*(T)X ] {T)dT - E ^ 2c k Xf(r)X k (r)dT 

Choosing now r € argmin[ s t j , we get 



>X){s) - E J 2c k X 2 k (T)X k ( T )dT 



By summation over all nodes j with \j\ — n we have 

E ™ ■ > E X '( s ) - / E 2c fc Jff (r)X fc (r)dr. 

|j| = n \j\ =n \k\—n-\-l 

Finally, we apply for m = n — 1 , n the following integral form of © to get the 
first part of the thesis. (Even if n = and m = — 1 this is true, trivially.) 

£rn(t) - £ m {s) = - f E 2c j X K T ) X j( T ) dT 

J S I ■ I II 

\j\—m-\-l 

We turn to the second part. Sum jj on every j with \j\ = n to get 

POO 

E 7; = E X * (°) + 7 2 E cjA-^da, 
by (fT4|) the integral term is bounded above by £ n _i(0), so 

E 7i< E ^ 2 (o)+^-i(o) = E^ 2 (°) = f ( )- 

n |fc|>n feE J 
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Finally, the third part. Let r > 0. By computing the time derivative of 
X)fceTi X% which is analogous to dropping the border term and integrating 

|fe|<n 

on [0, r], we have, 



x k(r) 2 < x k(0) 2 +2 [ 2c J X j X?du< 



7j 



k£Tj feGT. 
\k\<n |fc|<n 



Now, let n —> oo to conclude. □ 

The following statement will be used in the proof of Lemma 15.31 
Lemma 5.2. For every h > and A > the following inequality holds: 



I j e - x ^drds>^(l-e-^). 



Proof. 



o Jo 



< ] drds > / / e - A(s - r) dr ds = A fi , Mr). □ 



Lemma 5.3. Assume that a > a, where 2 2a = — $Oj is the constant 
number of children for every node. Let X be a positive I 2 solution of (|17p . Let 
(S n )n>o be a sequence of positive numbers such that J^ Tl S n and §- 2 2~'' a ~ a ' n 
are both finite. Then there exists a sequence of positive numbers (h n ) n >o such 
that h n < oo and for all n > for all t > 

£ n {t + h n ) ~ £ n - x {t) < 5 n . (18) 

In particular, for every M > 0, 

/ CO \ OO 

£\ Y\ K ) < £ M -i(0) + V S n . (19) 



The sequence 



n=M ' ra=M 



satisfies and (fH))) . 



Proof. Fix n > and positive real numbers i, /i„. For all j of generation 
n, let rrij :— inf,. e [ t l t+h n ] X 3 (r). We claim that if h n is defined by (|20[) . then 
S|j|=n m j — ^ni which together with Lemma I5TT1 completes the proof of (TT5T) . 

We prove the claim by contradiction: suppose that ;i_ „ rrij > S n . We will 
find a contradiction in the estimates on £ (0). By Proposition 13.41 

f (0) > 2 / ^ ^ c fc X fc (t + s)X?(t + s)rfs 

"' |j|=nfeeOj 

We have a lower bound for Xj , namely mj , but we need one also for X^ . 
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For all j € J, let Tj := max(7j,£ (0)iV* '"''). From Lemma 15.11 we have 
Ewi= n 7.j 5; ^(0) and hence Eh'^n^j — 2£(0); by the same lemma, for all 
i E Tj we have Xf < < Tj uniformly in time, so for all k E Oj, 

X k = C k Xj - ^ ClXiX k > CklTL* - XjX k 

where A,- = N*2 na+2a ,/f~. This eh 



X k (t + s) > c k m) J e- x ^ s -^dr 

We can write 

£(0) > 2 < f z- X][s - r) drds ]T c\ 
and by lemma [5T21 we have 



\j\=n J fceOj 



c 2 , 



Let us focus on the exponential. We substitute (|2U)) and make use of the in- 
equality Tj > £(0)N- n = £(0)2- 25 ™, 

^h n _ +2a ^ V2£(0)^ g(0) 2 

By the hypothesis that E|i|= n TO j ^ ^™ an< ^ Lemma lSTTl we know that <5„ < £(0) 
we get 1 — e~ A j'W 2 > I, We obtain 

Now we can use Cauchy-Schwarz and the AM-QM inequalities to get 

m) (E|,-|=„"*?) (E b |=„»^ 



E 



> =^=- > 



b1= 

again by the hypothesis that E|j|=n TO j > ^» an< ^ thanks to E|j|=n ^ < 2£ (0), 



A /2f(0)2 an 



so that the right-hand side of ([2"Tj) becomes larger than £ (0), which is impossible. 

We turn to the second part. Let M > and define the following sequence 
(t«)n>M-i by tM-i = and f n = t„_x + /i„. By (JTSJ with t = we get 

£n(t n ) — £ n -\(t n -\) < 6 n . 
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We sum for n from M to JV, yielding 

N 



n=M 

which, due to monotonicity of En, yields 

(oo \ N 

n=M ' n=M 

Now we let JV go to infinity to get the thesis. □ 

Remark 3. It is easy to prove this result also if relaxing the condition on the 
number of children from constant number to 1 < §Oj < JV*. One has to change 
slightly the definition of h n , which becomes 



£(0) 3 / 2 



S 2 , 



2-(a-a)ri+2a+3/2 



Theorem 5.4. Assume that a > a, where 2 2a = JV* = jJOj is the constant 
number of children for every node. Then for every e > and rj > there exists 
some T > such that for all positive I 2 solution of (jTTJ) with initial energy 
£(0) < i] one has S(T) < e. In particular 

lim £(t) = 

t— >oc 

i.e. there is anomalous dissipation. 

Proof. Given e > let us take a sequence of positive numbers (<5„) n >o such that 

oo oo _^ 

ES n = e and > — ; .„ < +oo 

n—o n—o n 

This is possible, for example, taking 5 n = e(l - 2- < - a - & ^ 3 )2- ( ~ a - & '> n / 3 . Now 
Lemma applies, so by the definition of h n given in (|20p 

h„ < ~, — rr — — and h n < ^~Jl — _ j> 

™ - 2 (a-a)«,52 ™ - (l - 2-( Q ^ Q )/ 3 ) 3 

n n=0 v ' 

Take M — in (fTO)) and by monotonicity of energy £(T) < e. □ 

We are finally able to prove Theorem 12.11 which is a consequence of Theo- 
rem 15.41 with a rescaling argument based on the fact that the non-linearity is 
homogeneous of degree two. 

Proof of Theorem \2.1l By Theorem 15.41 for every < p < 1 there exists r > 
depending only on p and £(0), such that £ (r) < p 2 £(0). We will apply this 
bound to many different solutions, all of which have energy at time zero not 
above £ (0). 
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Let $ = 1/p > 1. We can define the sequence 

( ? Qn _ 1 \ 
+ — — —t\. 



n > 1 



It is immediate to verify that all of these satisfy the system of equations (JTTJ) , 
but with possibly different initial conditions. We have 

E(^ n) (o)) 2 =^E( x i n_1) w) 2 - 

Recalling the definition of r, the above equation allows to prove by induction 
on n that for all n > one has £j eJ (X- n) (0)) 2 < 8(0). For all n > 0, let 

*» = 

Then by the definition of iW, we have proved £(t n ) 2 < i)- 2n £ (0). Since i? > 1, 
t„ f oo, hence given t > there is n such that t n < t < t n+ \. That means we 
have by monotonicity 

£{t) < tf- 2n £{0) and < \ 

finally, by definition t n+ i < \\ t — ^jz^,, so for C = £(0)( j-^) 2 we get 

£(t)<$- 2n £ (0)<-^<§ □ 

l n+l l 

5.1 Self-similar solutions 

We devote the end of this section to prove the existence of self-similar solutions. 
We call self-similar any solution X of system (fTT|) of the form Xj(t) = a,j(p(t), 
for all j and all t > 0. By substituting this formula inside p7|) it is easy to show 
that any such solution must be of the form 

for some to < 0- The condition on the coefficients aj is much more complicated 
a = 



< 



a? = ^ CkCLjdk, Vj e J 



fceo 



so we base instead our argument upon |2J , where it is proven existence and some 
kind of uniqueness of self-similar solution. We obtain the following statement. 

Proposition 5.5. Given to < there exists at least one self-similar positive I 2 
solution of (|17p with qq > 0. 
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Proof. We use Theorem 10 in [2] which, translated in the notation of this paper, 
states that there exists a unique sequence of non-negative real numbers (b n ) n >o 
such that 60 > and Y n := is a positive I 2 solution of the unforced inviscid 
classic dyadic ((5]) . Thanks to Proposition 14.11 this solution may be lifted to a 
solution of the inviscid tree dyadic with the required features. □ 

Remark 4. For the tree dyadic model self-similar solutions are many. In the 
standard dyadic case studied in [2] it is shown that given to < and no > 1 there 
is only one I 2 self-similar solution such that no is the index of the first non-zero 
coefficient. If no > 1, this solution can be lifted on the tree to a self-similar 
solution which is zero on the first n$ — 1 generations. We can then define a new 
self-similar solution which is equal to this one on one of the subtrees starting 
at generation no and zero everywhere else. Finally, we can combine many of 
these solutions, even with different no, as long as to is the same for all and their 
subtrees do not overlap. 

6 Stationary solutions 

In this section we will study the stationary solutions for both the classic dyadic 
model (|16[) and the tree dyadic one (TT5j) . We will in particular restrict ourselves 
to study positive I 2 solutions which are time independent. Proposition 14 . 1 1 allows 
us to link the two models, in that for any solution of the classic dyadic model 
one can build a solution of the tree dyadic model. Thus is it enough to prove 
existence for the classic dyadic and uniqueness for the tree dyadic. 

One purpose of this section is to prove the existence and uniqueness of sta- 
tionary solution on the tree dyadic model and extend existence and uniqueness 
results given in [3] and [5] for the dyadic model. In [6J it is proven that the 
inviscid dyadic model with ft = I has a unique stationary solution, while in the 
companion paper [7] it is proven that such a solution is a global attractor. The 
viscous dyadic model is studied in [5], where it is proven that for ft <= (§, §] the 
stationary solution is unique and is a global attractor. In [4] it is proven that 
for the viscous case it is possible, dropping the Y n > condition, to explicitly 
provide examples of non-uniqueness of the stationary solution. In this paper 
we prove the existence and uniqueness of stationary solutions in I for every 
positive value of the ft and 7 parameters both in viscous and inviscid dyadic 
models. This will provide a corresponding result of existence and uniqueness 
for a > a and 7 > in the tree dyadic model. Furthermore in the inviscid case 
we will explicitly provide those solutions (Proposition 16 . 1} . while in the viscous 
case we'll prove that the stationary solutions are regular if and only if is 
big enough, > 2 2a ~ 3 ~ f or the forcing term / is small. For f — the unique 
(non-negative) stationary solution is trivially the null one, so in this section we 
assume / > 0. 

6.1 Stationary solutions in the inviscid case: existence. 

In the inviscid case, the differential equation is very simple, so it is easy to find 
stationary solutions in the class of exponential functions. One immediately finds 
the following result. 
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Proposition 6.1. Consider the tree dyadic model (|15p and the classic dyadic 
model (|16p . both inviscid (v = 0). Let 2 2dl = N* = ttOj be constant for all j £ J. 
Then: 

1. the sequence of constant functions Y n (t) := / • 2~^(™ +1 ) is a positive I 2 
solution of the system (|16p . 

2. the family of constant functions Xj{t) :— f ■ 2 _ ^3^( 2a +") for j £ J is a 
positive componentwise solution of system (1 1 5 P it is also an I 2 solution 
iff a > a; 

Proof. A direct computation shows that X and Y are componentwise solutions. 
To show that Y is I 2 observe that, since /3 > 0, ||Y|| <oo. To check whether X 
is I 2 compute the energy by generations; we have for n > 0, 

|j|=« 

with C not depending on n. Hence X is ^ 2 if and only if a — a > 0. □ 

6.2 Stationary solutions in the viscous case: existence. 

In the viscous case, the recurrence relation coming from the definition of station- 
ary solution is more complex, and has no solutions in the class of exponential 
functions. Anyway, by careful control of the recurrence behavior, we are able to 
prove that a stationary solution exists, and also to distinguish if it is conservative 
or has anomalous dissipation. 

Definition 4. We say that a stationary positive I 2 solution X is regular if for 
all h € K 

Y^[2 hm Xj} 2 < oo (22) 
jeJ 

Theorem 6.2. There exists a stationary positive I 2 solution of the classic dyadic 
model (I16p when v > 0. 

Theorem 6.3. Consider any stationary positive I 2 solution of the classic dyadic 
model (|16p with v > 0. 

1. 7/37 > 2/3 then it is regular and conservative. 

2. If 37 < 2/3 then there exists some C > such that if f > C the stationary 
solution is not regular and there is anomalous dissipation. 

Before we go into the proofs of these theorems, let us introduce a useful 
change of variables, that will come handy in both proofs. If Y is a stationary 
solution of (TIT))) then, for every n > 0, we have 

-vT< n Y n + 2^Y 2 _ X - 2 f3n +?Y n Y n+1 = 

This equation can be made into a recurrence, and the change of variables that 
best simplifies its form is 

Z n := ! y- 1 2K"+ 2 )y„. (23) 
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Since the stationary solution in the inviscid case decreases like 2~ a"", the expo- 
nent's rate is in some sense expected. The system of differential equations 
for Z becomes 

f Z_ 1 = ^ 1 2l/=:.g 

Proof of theorem \6.2[ Let us consider the change of variable ([23)) , we have to 
show that the system has a positive solution for which Y is I 2 . System (flM)) 
gives a recursion which, given Z_i = g and Zq allows to construct the sequence 
(Z n ) n >_i in a unique way. Any such sequence will give a stationary componen- 
twise solution. What we want to prove is that there is some value of Zq such 
that this turn out to be a positive I 2 solution. Let we exploit the dependence 
from Zq by defining a sequence of real functions 

Z_i(a) = g 

Z (a) = a (25) 

z «+i( a ) = %Sr- 2(7 ~ §/3) "< n >° 

Now we construct a descending sequence of open real intervals (I n ) n >o such 
that (0,oo) = Jo D Ii D 1% D . . . and such that Z n is continuous and bijective 
from I n to (0, 00), with Z n strictly increasing for even n and strictly decreasing 
for odd n. 

Let Iq — (0,oo). Z{j{a) is monotone increasing, continuous and bijective 
from Jo to (0, 00). 

By (|2"5|) we have that Z\(a) — g/a 2 — 2^ 7 ~3^) is monotone decreasing, con- 
tinuous and bijective from Jo to (— 2 ( - 7_ 3< 9 ) ; 00) so there exists a limited interval 
(61, ci) := I\ C Jo such that Z\{a) is monotone decreasing, continuous and 
bijective from 1\ to (0,oo). 

Now suppose we already proved for m < n that Z m {a) is continuous and 
bijective from I m to (0, 00), with Z m strictly increasing for even m and strictly 
decreasing for odd m. 

Suppose that n is odd (resp. even). Then by (|23|) Z n+ i(a) is monotone in- 
creasing (resp. decreasing), continuous and bijective from J„ to (— 2( 7- 3 /3 '™, 00) 
so there exists an interval (6 n +ij c n+i) '■= J1+1 C I n such that Z n+ i(a) is mono- 
tone increasing (resp. decreasing), continuous and bijective from J„+i to (0,oo). 

Observe moreover that the borders of these intervals are not definitively 
constant, since for all n, 6 n +2 7^ b n and c n +2 7^ c„. Hence if we define b — lim n b n 
and c = lim„c„, it is clear that for all n, b n < b < c < c n , that is the closed 
interval (possibly degenerate) [b, c] is contained in every J„ . 

Now we choose any a £ [b, c] and we know that the sequence Z n (a) is strictly 
positive. We are left to prove that it is also I 2 . To this end let Y n be any 
stationary, positive componentwise solution. Let £ n = X)fc=o ^fe m ana l°gy 
with the definition for the tree model. We compute the derivative 

= j t £ n {t) = -v^2 l k Y 2 + f 2 Y - k n+1 Y 2 Y n+1 , 

k<n 

hence, since l k >l,£ n < J2k<n lkY k ^ ^ 1 P Y o for all n. □ 



Vn > 



(24) 
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Proof of theorem \6.3[ Let us consider again system (|24[) and let \i := 7 — If 
> the corrective term goes to infinity, while if fj, < it goes to zero, so we 
expect two different behaviors in the two cases. We'll show that in the first case 
Z n goes to zero super-exponentially for n — > 00, while in the second one Z n \. z 
and z > if g is large enough. 
Case n := 7 - |/J > 0. From ([HJ we get 

nun ^2 ^2 v ^2 17 

A 6 n — /j n _ x £ n — A n A n+ i 

Sum over n to get 

J2 ^ k Z 2 k = g 2 Z - Z 2 n Z n+1 (26) 

k<n 

Since [i > 0, by positivity of Z, we have 

lim Z n = 0. (27) 

n—too 

From (T23]) and Z n+ i > we get Z n < Z^_ x and since by (gTJ) Z fi =: A < 1 for 
some n, by iterating the above equation we get for all m > 

Zn+m < A 2 

that is to say that Z n goes to zero for n going to infinity like the exponential of 
an exponential, so for every s > we have 

J2 {2 sn Z n f < +00 and J2 ( 2SnY nf < +00. 

n n 

It is now clear that lim„ k n+ iY^Y n+ i = 0, so Y is conservative by Definition [3] 
Case /1 := 7 — §/? < 0. The first step is to prove that Z n is non-increasing 
in n. Suppose by contradiction that for some n we have z Zn 1 = A > f , then we 

claim that |^ > A 4 > f and hence by induction Z Z "+ 2 " > A 4 ™ . By (TJU) for 
all k > 

r7 . Z k-l _ Zk-l ry 

^k+1 < — — — &k-l 
^k ^k 

This can be used iteratively together with the claim to show that 
7 2 

ry _ * J n+2m-l r,a(n+2m) 

^n+2m+l — — y 6 

^n+2m 

Z n +2m-l Z n+2m -3 Z n -\ _ ^{n+lm) 

ry ry ' ' ' ?l 1 

^n+2m ^n+2ra-2 &n 

< Z n ^\- im - 2^(™ +2m ) 

so we get a contradiction because Z„ + 2m+i < for some m. 
We prove the claim. Let x — z $ n = z . Observe that 

V"2 yfin 

Z n+1 = ^-2^ = —(l-x) (28) 
We divide by Z n (and we notice that x < 1), 

^=\- 2 {l-x) (29) 
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Now 

Zn+2 = _ 2^(«+i) > Z 



so dividing by Z n +\ and substituting (f2U|) and ([25]) . we get 

Z n+ i X n+ i 1 - x 1 - x 

Since A > 1 > x > 0, it is now clear that > A 4 . So we have proven the 

claim and showed that {Z n } n >o is non-increasing in n. 

The last step is to show that for g large enough Z n I z > 0. By rearrang- 
ing (|2"6"]l and recalling what we proved above, 

n n 

Z 3 > Z 2 Z n+1 = g 2 Z Q - J2 ^ k Zl > g 2 Z a - gZ £ 2» k > gZ (g - ^—^j 

so if g > j-j^g then Z„ converges to a strictly positive constant z. 
To prove anomalous dissipation we compute the limit 

lim k n+1 Y 2 Y n+1 = lim 2^+^ 3 2~^- 7 ^ 3 ^Z n+1 = 2~ 4 ^Vz 3 > 

n— >oo n— >oo 

So by Definition [3] there is anomalous dissipation. □ 

6.3 Stationary solutions in the inviscid and viscous case: 
uniqueness 

We prove uniqueness in the class of stationary positive I 2 solutions for the tree 
dyadic model. The result also holds for the classic dyadic, because it is a 
particular case of the former, or by virtue of the lifting Proposition 14. II 

Theorem 6.4. Consider the tree dyadic model ([7]) and assume that a > a, 
where 2 2a = N* = §Oj is the constant number of children for every node. Then 
there exists a unique stationary positive I 2 solution. 

Proof. Existence is a consequence of Proposition ^, ll in the inviscid case (v = 0) 
and Proposition 14. II and Theorem 16.21 in the viscous case. 

To prove uniqueness we apply a change of variables similar to (|23p 



(2+Mla 

7., : 2 ' A,.. VjeJ (30) 



Then from ([7]) we have 



jZj = -vV^Zj + 2% a ^Z 2 - ^ aljl ZjZ k (31) 

so if X is a stationary solution, Z must satisfy 

Z- = f ■ 2 Q / 3 

Y j Z k = ^--v2^~^\. (32) 
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Moreover observe that the condition X <E I 2 is equivalent to 

]T(2-§l^-) 2 <oo (33) 

Assume by contradiction that there are two different stationary solutions 
of (p?2"j) which we denote by W = {Wj}j^j and Z = {Zj}j e j. Let n be the 
smallest integer such that there exist j\ G J with |ji| = n and Wj 1 7^ 
Without loss of generality we can take ^r 1 - =: A > 1. 

Let jo = fco — ji and fci = j%. Extend these to two sequences of indices 
(j m )m>Q and (fc m ) m >o with j rn e O jm _ 1 and k rn € Ofe m _ 1 , picking alternatively 
among those that maximize or minimize Wj m and . 

More precisely for m > 2 choose j m € Oj m _ 1 and fe m € O km _ 1 in such a way 
that if to is even 

Wj m = min{Wj : i € Oj m _ 1 } Zfc m = max{^ : i £ Cfc m _ 1 } 

and if m is odd 

W Jm = max{Wi : i € Oj m _J Z fcm = min{Z 4 : i e Ofc m _J 

The idea supporting the definition of these sequences is to choose the indices so 
that 

W h < Z kl , W h > Z k2 , W h <Z k3 , 

We will now prove that, with our construction, those inequalities hold and, 
moreover, the ratio between W m and Z m grows according to 



> Jm - 1 • ' > A Vm > 2 even (34) 

Jm-2 



w jm - z km _, w: 



^-vfc'<I >A V 777 >30dd. (35) 

We prove inequalities (13"4"|) and (|33|) by induction on m > 2. First note that for 
777 = and 777 = 1, 

§^ = 1 and ^ = A (36) 
•Tjo Z fei 

Now we proceed by induction. Let 777 > 2 even. By the definition of j m , k m and 

by fl32J) we get 



W jm = min Wi <N- X W i = 



Z*. = max Zi > N^ 1 V Z t = N~ 



w 2 



(37) 



Z 2 



v 2(i-%ot){n+m-2) 



(38) 

By (I36p when m = 2 or by inductive hypothesis (fM)) and ([33]) when 777 > 4, 

/^L i = fL^^ i> /a m = 2 

3^-1 / Wj m _! W? Z km _ t ~ \ ( A 2 m - 4 )2 A 2™- 3 = A 2- 2 m > 4 
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so in particular the ratio is above 1 and, since for every a > b > c > we have 

j,2(7-fa)(n+m-2) 



§— - > f , for m > 2 even 

o— c — b 5 — 




Wj m w ^-2 2 (^f «)("+ m - 2 ) 

This concludes the inductive step for m even; for m odd the reasoning is analo- 
gous. We now want to use inequalities and ([33]) to get a contradiction. We 
will consider separately the cases v > and z/ = 0. 
Case v > 0. Let m be even; by (|34|) 



■Z 2 2 VF 2 

fcm-2 \ 2 

applying (f3"2")l to Wj m _ 1 we have 

VF 2 



Jm-2 



so from (I38p . putting everything together, we get 



k ^-2 _ J/2 (7-§a)(n+™-2) 



> iY-V2^-i Q )("+ m - 2 )(A 2 "^ 2 - 1) 



For m even going to infinity we have obviously that Z^ m grows as the exponential 
of an exponential, which is in contradiction with (1331) . 

Case v = 0. If v = we already know one explicit stationary solution, by Propo- 
sition UTTl namely Xj = f2~ ]j ^~ ( 2a + Q ). By the usual change of variables (|50|) 
Vj = /23( a_a (ljl -1 )) is a solution of (|3"2"j) satisfying the regularity condition ([3"3"| . 
Without loss of generality we can suppose that Wj — Vj or Zj = Vj . In the first 
case, for m even 

z k m > W jm \ 2m ~ 2 = /2§( Q - 5 (»+ m - 2 )A 2m ~ 2 
in the second case for m odd 

In both cases the right-hand side grows super-exponentially as m — > oo and this 
is in contradiction with (1331) . □ 



Proof of Theorem \2.cft Existence and uniqueness are given by Theorem 16.41 

If v = the solution is identified by Proposition HTTJ If v > 0, by uniqueness, 
the solution is the lift of the stationary solution of the classic dyadic with j3 = a— 
a, as per Proposition ^. II Then the two regimes are proven in Theorem l6.3l □ 
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